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Latest Results from Heavy Quark Simulations 

John H. Sloan 

SCRI, Florida State University, Tallahassee, Fl 32306-4052, USA 

I review the status of bh and cc calculations, numerical and analytic. I discuss the extraction of Os and quark 
masses from spectrum calculations. I compare the NRQCD and Improved Heavy Wilson formulations of heavy 
quarks, and discuss recent calculations using a Heavy Staggered formulation. 



1. Introduction 

I have taken my charge in this review to be 
two-fold: to review this year's results from heavy- 
heavy meson calculations and to compare the 
two main formalisms for treating heavy quarks, 
NRQCD [| and the mass dependent improve- 
ment program of the Fermilab group (which 
I call Improved Heavy Wilson (IHW) fermions). 
I will try to minimize my overlap with previ- 
ous reviews |^-|^. The messages that I hope 
people take from the theory part of this talk 
are that IHW fermions are (improved) Wilson 
fermions (with an on-shell renormalization con- 
ditions), that NRQCD and IHW differ only by 
negligible irrelevant operators at large am, and 
that NRQCD actions have an extra symmetry 
(the heavy quark symmetry) which makes them 
much less messy to deal with but causes them to 
blow up at small am. 

From the new results, my message is that we 
probably understand what we're doing in onium 
systems; discretization errors, quenching effects, 
spin splittings and perturbation theory all seem 
to behave as expected. If this is true, then we are 
probably within one to three years of being able 
present an extensively cross-checked determina- 
tion of a^{Mz), with uncertainty at the 1-2% 
level and all ingredients calculated in several dif- 
ferent ways, to the rest of the particle physics 
community. Another message is that tadpole- 
improved coefficients Q are necessary for pertur- 
bative improvement to work well. Finally, there 
is a discrepancy between typical lattice spacings 
determined from bb physics and those from light 
spectroscopy (both quenched and dynamical); we 



should be aware of this as a potential source of 
systematic error in light hadron calculations. 

The reader should assume that any result re- 
viewed here is preliminary, unless it appears else- 
where. 



2. IHW and NRQCD: Compare, Contrast 

2.1. A Toy Model 

The most important feature of heavy quark 
(am^l) calculations on the lattice is that, for 
a lattice 0(4) invariant action, boost invariance 
of on-shell correlation functions is strongly bro- 
ken, even for soft boosts of particles with (spa- 
tial) momentum much smaller than the cutoff. 
This is because all components of the on-shell 
four-momentum cannot be made small simulta- 
neously; the energy will always be of order the 
cutoff or higher. To restore low-momentum boost 
invariance without taking the cutoff much larger 
than the mass, (lattice) boost invariance must 
be broken in the bare action; lattice 0(3) is 
the appropriate symmetry group for heavy quark 
bare actions. NRQCD implements the breaking 
0(4)^0(3) directly; NRQCD is a non-relativistic 
effective action. In the IHW formalism the break- 
ing is more subtle; the IHW formalism can be 
thought of as just the normal improvement pro- 
gram for Wilson fermions, except the renormal- 
ization conditions are imposed with the heavy 
quark 0K-s/ien|Q. The Lorentz breaking arises 
because the four-momentum of a heavy quark 
is very anisotropic; a 'preferred frame' is chosen 
which is just the rest frame of the heavy quark. 

To illustrate this, consider the dispersion rela- 
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tion of an on-shell (free) naive fermion; 

(am)^ + ^sin^(ap^) = 0. (1) 

A" 

If am <C 1, all of the components of the dimen- 
sionless momentum can be chosen small, so that 
sm can be replaced with its expansion; 

(am)2 + (ap)2-i^(ap^)4 + ... = 0, (2) 

where {ap)'^=^^{apf^)^ . The first two terms in 
(§) represent the continuum dispersion relation, 
while the third and later terms are cutoff effects. 
Since these are much less than one, (||) is approx- 
imately Lorentz invariant at low momenta. 

When am is large, however, this fails. At least 
one component (the energy) of the 4-momentum 
in (|^) must be large, so as to cancel the mass 
term, and the sin of this component cannot be 
expanded. This leads to the dispersion relation 

{amf+sm^{apo)+{apf~^"^{ap^)^+. . . = 0, (3) 

i 

where p is the spatial momentum. Even at low 
momenta, (||) is not Lorentz invariant, and so 
the bare action that leads to it cannot be used 
directly to calculate S-matrix amplitudes. Dis- 
cretization errors in the temporal and spatial di- 
rections are different; to obtain the correct behav- 
ior near the mass-shell the temporal and spatial 
terms in the action must be adjusted separately. 
This results in an action which only has lattice 
0(3) symmetry, but which has 0(4) invariant cor- 
relation functions at low momentum. 

The dispersion relation (^ has behavior very 
similar that of IHW fermions ; the static 
mass. Ml, depends logarithmically on am when 
am ^ 1, while the kinetic mass, M2, goes lin- 
early with am. A Lorentz invariant theory has 
Mi=M2=M3—. . ., where M3 is the mass appear- 
ing in the p** term in the relativistic low momen- 
tum dispersion relation. Mi=M2 can be recov- 
ered by rescaling the spatial derivatives relative 
to the temporal derivatives; this is equivalent to 
the 'asymmetric k' prescription in ||^. M3 must 
be adjusted by adding a p^ improvement term to 
the action; note that this is a lattice 0(3), rather 



than lattice 0(4), invariant operator. For Wilson 
fermions, the clover term, 'ijjaf^^F'^'^ilj, is the low- 
est order improvement operator needed in the ac- 
tion ||^. When am 3> 1, the asymmetric physical 
renormalization point requires that a ■ F he bro- 
ken up into a sum of electric and magnetic 0(3) 
invariant representations, with separately tuned 
coefficients for the electric and magnetic terms. 
I will call this the "asymmetric clover" level of 
improvement for the IHW formulation. 

The Fermilab group has pointed out that 
only a modified Lorentz invariance is needed 
when studying processes involving valence heavy 
quarks; rather than using an asymmetric k to 
tune M1—M2, one can shift the zero of energy 
in the continuum theory by Ai2 — Mi for each 
quark without changing the physics. If this is 
true, then k can be kept symmetric as long as Mi 
is ignored. This is typically done both in IHW 
and NRQCD. Asymmetric improvement terms 
are still needed at higher orders; for example 
an asymmetric clover coefficient is needed to get 
both spin-spin and spin-orbit interactions right. 

2.2. Similarities and Differences 

IHW and NRQCD are very similar formalisms 
for simulating heavy quarks. Both break lattice 
0(4) invariance to lattice 0(3) in the bare ac- 
tion to recover boost invariant low momentum 
correlation functions. Both rely heavily on an 
improved action program to eliminate discretiza- 
tion errors, and the systematic errors in both 
can be expanded in Aqcd/Mq (where Aqcd ~ 
SOOMeV). Although the IHW is formulated in 
terms of Dirac spinors, it can be thought of in 
the Foldy-Wouthuysen-Tani(FWT) transformed 
basis of NRQCD. 

There is an intrinsic difference between the two 
formalisms, however. In the FWT basis, NRQCD 
throws away the coupling between large and small 
components of the 4-component spinors, resulting 
in a block diagonal action involving (2 sets of) 2- 
component Pauli spinors. This imposes an addi- 
tional symmetry on the action (the heavy quark 
symmetry), and causes the NRQCD action to be 
non-renormalizable when am —f 0. The advan- 
tage of this procedure is that m can be removed 
as a scale before discretization. The IHW action 
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does not eliminate the couplings between large 
and small components and so is renormalizable 
(going over into an improved Wilson action as 
am ^ 0) , at a cost of having to treat discretiza- 
tion errors to all orders in am. 

To illustrate the difference, start with the 
free continuum Dirac action (note that neither 
NRQCD nor IHW is actually derived as below, 
but the basic principles are the same); 



dt 



icr ■ p 



icT ■ p —dt + m 



V'- 



(4) 



A continuum FWT transformation can now be 
performed to some order in v'^—{'p/m)'^, resulting 
in (for the lowest order transformation) 



- / dt+m + ^ 



dt + m^^ 



(5) 



where . . . represents terms which lead to higher 
order corrections to the dispersion relation, such 
as p^/m"^. If we discretize (||) directly, we obtain 



a ^ sinh(a£^) + m + ^ + . . . 
— a^^ sinh(ai?) + m + ^ 



(6) 



This is the IHW approach; as m gets large, 
sinh(a-B) and hence aE must also become large, 
leading to large discretization errors. In the 
NRQCD approach, the higher order terms (. . .) 
in (H) are set explicitly to zero. This makes (0) 
block-diagonal; the transformation {[) 'ijje"^'^°^, 
ip e~™''°*tp can now be performed on (^, re- 
sulting in an action 
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(7) 



This is the NRQCD approach; ^ is 0{mv^), 
so there are no problems discretizing (|^). When 
am — > 0, the symmetry imposed by the zero en- 
tries in (|^) is no longer valid; the theory cannot 
be renormalized. This method for removing the 
zero of energy doesn't work for IHW; 70 doesn't 
commute with the off-diagonal components of (||) . 

2.3. Pros and Cons 

In IHW, discretization errors must be treated 
to all orders in am. In NRQCD, m has been 



removed as a dynamical energy scale, and dis- 
cretization errors can be expanded in ap. This 
means that: 

• In IHW, the tree-level coefficients of terms 
in the action are complicated, non-intuitive 
functions of am. In NRQCD, these coeffi- 
cients are simple functions, for which non- 
relativistic intuition can be easily applied. 

• NRQCD is a non-renormalizable theory; 
the coefficients in the NRQCD action blow 
up when (see Fig. |^). An IHW ac- 
tion is renormalizable; as am — > it goes 
smoothly over to a perturbatively improved 
Wilson action. 

• Because NRQCD blows up at small am, 
charm can only be simulated at /3^5.85. 
IHW simulations can be performed at any 
/3; brute force can be used to reduce dis- 
cretization errors and there is a much larger 
range of (3 for which charm and bottom 
quarks can be simultaneously studied. 

• Because of the complicated coefficients, 
the IHW action is harder to derive than 
NRQCD. To date, only the terms which 
contribute at leading order in v'^ to the spin- 
averaged and spin dependent spectra have 
been published for the tree-level IHW ac- 
tion. The relativistic corrections to these 
terms have been published, also at tree 
level, for NRQCD. 

• IHW fermions cost more CPU time than 
NRQCD; IHW propagators solve an ellip- 
tic partial differential equation(PDE), while 
NRQCD solves a parabolic PDE. I sus- 
pect the IHW inverter can be accelerated. 
When am is large, the IHW PDE is close 
to parabolic; there should be a way to use 
this. 

2.4. Levels of Improvement 

I find it easiest to use the NRQCD action de- 
rived in ^ to catalogue improvement terms in the 
NRQCD or IHW actions and to understand their 
effects. The 0(mw'*) NRQCD action is given by 

S = ij{At + H)t^. 
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where A*^^^ is the discrete laplacian, and 
A(*)=^. A|. For HH mesons, the d terms all 
contribute 0{mv'^) to the energy; the kinetic en- 
ergy term is, of course, O(mv^). Further improve- 
ment operators are suppressed by more powers of 
~ Kqcd/M, which is about .1 for bb and .3 
for cc. The counting for improvement operators is 
different in HL mesons; see ||] for details. All op- 
erators are tadpole- improved; Ci=l at tree- level. 

The terms in (|^) are easily identified from low 
energy electrodynamics; Ci multiplies the leading 
relativistic correction to the dispersion relation, 
C2 the Darwin term, C3 the spin-orbit coupling 
(causing fine-structure), and the C4 term gives 
rise to the spin-spin interaction (causing hyper- 
fine splittings). The C5 and Cg terms correct 
the leading discretization errors in the spatial and 
temporal directions, respectively. 

As more accuracy and physical information is 
required from a heavy quark calculation, more 
improvement terms can be turned on in the ac- 
tion. For every NRQCD action, there is a cor- 
responding IHW action 1^ in which similar im- 
provement terms have been turned on. At large 
am, the effects of improvement terms which have 
not yet been turned on should be the same size 
in IHW and NRQCD. At small am, however, 
the IHW action recovers lattice 0(4) invariance, 
so the effects of lattice 0(4) breaking improve- 
ments are small. In NRQCD language, the Wil- 
son action gets the term corresponding to Ci right 
at small am, without needing improvement. In 
this regime, it is much harder to estimate the 
IHW systematic errors and compare to O(mw^) 
NRQCD, for example. El-Khadra has used po- 
tential models to estimate these errors in the Fer- 
milab group's calculations |^ , but it is also impor- 
tant to obtain simulation results from the IHW 



action equivalent to (||) , so that the two methods 
can be directly compared. 

The simplest action to calculate with is just 
0{rav^) NRQCD, in which all the Ci terms in 
(^ are turned off. This is a spin-symmetric ac- 
tion; spin-averaged splittings will have errors at 
the level, but spin dependent splittings will be 
zero. The equivalent IHW action is just the nor- 
mal Wilson action; spin splittings are non-zero 
but much too small. Note that the Ci and C2 
terms contribute to the kinetic mass at the 
level; if more precision is needed for A/2 then these 
improvements (and probably C5 and Cg) should 
be turned on. 

To study spin splittings, the C3 and C4 im- 
provements must be turned on. In IHW, this is 
equivalent to turning on the clover term; the mag- 
netic part corresponds to the C4 term, while the 
electric part is the sum of C2 and C3 improve- 
ments. The Sheikholeslami-Wohlert action|^ has 
the correct C4 improvement, but C3 and C2 are 
somewhat wrong. This is fixed by going to the 
asymmetric clover action; at tree-level, the C2, 
C3, and C4 terms are all correctly included. 

The most complicated action in use is the mv'^ 
NRQCD action in (^). Errors in spin averaged 
splittings and kinetic masses should be at the 

level, while spin splittings still have errors. 
The corresponding IHW action is not published; 
it is important that IHW calculations be done at 
this order to compare with NRQCD results. 

If one wants to reduce errors in the spin split- 
tings to 0(1;^), then one must go to the level 
of improvement in [Q. This is the first level 
at which one-loop results are needed for the Ci. 
Since as should be about the same size as , the 
aginv'^ errors arising from the use of tree-level Cj 
should be about the same size as mv^ improve- 
ment terms. (Note that I have used this to sim- 
plify the discussion of errors above, v'^ errors re- 
ally means v"^ and asv"^ errors; both should be 
the same size.) Until the one-loop corrections to 
C3 and C4 (or their IHW equivalents) are known, 
this level of improvement cannot be fully imple- 
mented. 

The results presented this year have differ- 
ent levels of improvement for different groups. 
The mw^ NRQCD action was used by the 



5 



UK(NR)QCD group|ig], and by Bodwin, Kim 
and Sinclair who have performed a very nice 
calculation of matrix elements involved in bb de- 
cay. The KT group and Wingate et a/.|l|] 
(referred to as KT and CDHW below) both use 
a Wilson action, which should have similar er- 
rors. The Fermilab group]^ uses the SW action; 
the NRQCD collaboration ^^^\ uses the mv 
NRQCD action. Both use tadpole improvement, 
which is essential. 

3. Perturbation Theory 

Perturbative calculations are an essential ingre- 
dient in both the NRQCD and IHW formalisms. 
To date, NRQCD and IHW simulations have been 
done using only tree-level coefhcicnts for the im- 
provement terms in the action. The next level 
of relativistic corrections, O(rnu^), however, are 
about the same size as the one-loop corrections 
to the current action, 0{asmv'^), so loop calcu- 
lations are needed to improve the actions beyond 
0{mv^). Besides the action, perturbation theory 
is needed for non-spectroscopic results. For exam- 
ple, a calculation of the pole or MS quark masses 
requires a perturbative calculation of the quark 
dispersion relation, while wave function and cur- 
rent renormalizations will be needed for quanti- 
ties such as fg or the leptonic width of onium. 

The Fermilab group has an ongoing program 
to perform perturbative calculations for the IHW 
formahsm[|jl|-|20|, while Morningstar has per- 
formed most of the recent perturbative calcula- 
tions in NRQCD [|l],|2|. Both calculate the dis- 
persion relation of the heavy quark; matching to 
the low energy continuum dispersion relation then 
fixes the Ci, C5, and Cq improvements, as well as 
the corrections needed to determine the pole or 
MS mass|17 . The main new results this year are 
the inclusion of the Lepage-Mackenzie q* in these 
perturbative calculations [^,| 19 1. 

The corrections to the action appear well be- 
haved in both formalisms; Morningstar's results 
for the one loop contribution to the C"s are shown 
in Fig. |l]. The one-loop correction is roughly 10% 
for am > .8; this is a good indication that pertur- 
bation theory is working. Notice that the correc- 
tions begin to grow rapidly below am < .6; this 
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Figure 1. Calculation by Morningstar [g2| of the 
1-loop perturbative correction to the coefhcicnts 
Ci, C5, and Cq in the NRQCD action as a func- 
tion of aM. Note that Morningstar's ci corre- 
sponds to a combination of Ci and Ce, while C2 
is C5 in d). 



is where NRQCD begins to "blow up." 

The Fermilab group has similar results for the 
perturbative correction to the static mass; their 
results for both SW and Wilson actions are shown 
in Fig. ^. Again, the result is very well-behaved. 
This calculation can be combined with simulation 
results to determine the quark's pole mass||]. 

4. Heavy Staggered Fermions 

This year, the KT group has introduced a third 
heavy fermion formalism; Heavy Staggered (HS) 
fermions ||2^ . They have performed simulations 
at both quenched (3=6.0 and at n/=2 staggered, 
/3=5.7, ams=. 01 11^, with bare valence masses 
ranging from amQ=.2 to amQ=.5. They used 
static, rather than kinetic, masses to match with 
experiment, but at these values of the am this is 
probably not a huge effect. As with most stag- 
gered fermion calculations, they had to deal with 
a daunting number of meson creation operators; 
they used 128 operators grouped into 36 irre- 
ducible representations! 

Results for the cc spectrum (averaged within ir- 
reducible representations) are shown in Fig. ||. As 
can be seen in the figure, the masses of the states 
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Figure 2. Perturbative correction to the static 
mass (Ml) for both Wilson (c=0) and SW (c=1.4) 
actions as a function of the kinetic mass (M2) j23|. 



Figure 4. Finite volume dependence of HS cc 
states. Note the large shift in the he state. Q 
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Figure 3. Heavy staggered cc spectrum. NG 
refers to flavor singlet pseudoscalar. Different 
points correspond to different irreducible lattice 
representations of the continuum statep4|. 



change in the IS-IP splitting. The smaller lat- 
tice is roughly the same physical size as most of 
those which have been used for cc simulations; an 
explicit check of volume dependence in cc spec- 
troscopy should probably be done using the IHW 
or NRQCD formulation. This is much less likely 
to be a problem in bb spectroscopy; the bb mesons 
are much smaller than cc. 

The KT group has extracted a value of 
a"^{Mz) from the cc IP-IS splitting; they ob- 
tain .111(4) (2) from their quenched results and 
an upper bound of .119(3)(3) at nf=2, where the 
first uncertainty is from the scale determination 
and n f correction and the second comes from fla- 
vor breaking. They have asked me to stress that, 
because of the strong finite volume dependence 
of the he mass, the nf=2 result should only be 
considered an upper bound. 



are dependent on the particular operator used. 
Furthermore, this operator dependence seems to 
go like the square of the heavy mass|^^. The 
most pronounced of these splittings is between 
the fiavor singlet and fiavor non-singlet versions 
of the Jjc', the KT group quotes separate results 
for the 'NG' and 'non-NG' pseudoscalars. 

A result that may be of concern to those us- 
ing other formulations for cc is shown in Fig. ^, 
which compares HS results at quenched (3=6.0 
on 16^x40 and 24^x40 lattices. The change of 
about 3% in the mass of the he results in a 20% 



5. Onium Spectroscopy 
5.1. bb 

Three groups have calculated more than one 
splitting in the bb spectrum recently. This year, 
the Fermilab|§,|] and NRQCD^,|| groups 
have increased their quenched statistics and the 
NRQCD group has results at n/ = 2. The results 
for spin averaged splittings are shown in Fig. |^. 
The Fermilab group determines the lattice spac- 
ing from the IP-IS sphtting; the NRQCD collab- 
oration uses both the IP-IS and the 2S-1S split- 
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ting. The zero of energy in Fig. || is adjusted so 
that the T(1S) energy matches experiment; the 
bare mass of the b quark is adjusted so that the 
kinetic mass, M2, of this state agrees with its ex- 
perimental value. 



GeV 
10.5' 



10.0-- 



9.5-- 



-i— -3S 



NRQCD (nf^2 
NRQCD [nj = 
UK(NR)QCD 
□ : Fermilab 



ID 



■f2S 



T 



VbiD) 



Figure 5. Spin averaged bb spectrum. A 
horizontal line indicates the experimental en- 
ergy; experimental values for spin-singlet P states 
are assumed to equal spin average of the cor- 
responding triplet states. Results are shown at 
/3=5.6 (n/=2 NRQCD), 6.0 (n/=0 NRQCD), 6.2 
(UK(NR)QCD), and 6.1 (Fermilab). 



The discretization and relativistic errors in 
Fig. H are expected to be about 5MeV for the 
NRQCD group, 50MeV for UK(NR)QCD[0, 
and somewhere between the two for Fermilab. 
This is because the NRQCD collaboration in- 
cludes all O(mw^) terms in the action, while 
UK(NR)QCD only works at 0{mv'^). The clover 
action used by Fermilab has an incorrect coeffi- 
cient for the p'^/(8rn^) term in the action, which 
becomes correct as am — > 0; the error in this co- 



efficient controls the size of O(tow^) errors in the 
Fermilab results. 

The NRQCD collaboration has spectrum re- 
sults at nf=2. These use an ensemble ob- 
tained from the HEMCGC collaboration, with 
two flavors of staggered fermions at am=. 01 and 
/?=5.6|2§. The CDHW results, which 1 will dis- 
cuss later, used this ensemble and a similar one 
(also from HEMCGC) at aTO=.025. The only 
other dynamical results available are from the KT 
group, who worked at P—5.7 rather than 5.6p5[. 
Both the CDHW and KT groups have results only 
for the fP-lS splitting. 

The HEMCGC collaboration has performed a 
tremendous public service by making their en- 
sembles publicly available; given the high cost of 
generating dynamical ensemble, it is important 
that as much analysis as possible be done on each 
one. Many other groups have been similarly gen- 
erous in providing their quenched ensembles to 
the NRQCD collaboration. 

The results show excellent agreement with ex- 
periment. The most important thing to focus on 
in Fig. ||is the relative size of the 2S-1S and f P-IS 
splittings; their ratio is plotted in Fig. ^ as a func- 
tion of . ft appears that there is a discrepancy 
in this ratio at n/=0, (2S-f S)/(f P-lS)=f .41(5), 
when compared with the experimental value of 
f .28, and that this discrepancy is removed when 
extrapolated to an Uf value between 2 and 4. 
Fig. H also shows consistent, but noisier, results 
at (3=5.7. This is an encouraging indication of 
scaling, but it should be noted that this splitting 
has been adjusted by a large («40MeV) pertur- 
bative correction which removes O(a^) errors in 
the gluonic propagator ||,0; The NRQCD col- 
laboration results may also show a quenching ef- 
fect for the 2P-fP splitting; this is not supported, 
however, by the UK(NR)QCD results. 

Fig. ^ exhibits exactly the behavior one expects 
from quenching errors; determining a^^ from the 
fP-fS splitting is roughly equivalent to match- 
ing the force in the quenched static potential at 
the Bohr radius to be the same as the force in 
full QCD, i.e. anj=o{ra)=afuii{ro) (this is the 
essence of the Fermilab prescription for correct- 
ing quenching errors in Ug determinations ) . 
Since the quenched /3-function runs faster than 
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Figure 6. Dependence on n/ of the ratio of the 
25'- 15 to IP - 15 splittings in bb. The filled cir- 
cles give values at quenched /? = 6.0 and (3—5.6, 
nf=2. The open circle is quenched /3 =5.7 for 
comparison. The horizontal dashed line gives the 
experimental value of 1.28 and solid line is an ex- 
trapolation to Uf — 3, drawn to guide the eye. All 
results are from the NRQCD collaboration 151. 
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Figure 7. hb spin splittings. The expected (non- 
quenching) systematic error in the NRQCD calcu- 
lation is shown; the Fermilab results should have 
similar errors. Symbols are the same as in Fig. |^. 



the full one, the potential a{r)/r will be weaker 
than it should be at short range. This in turn 
should shift quenched 5 states up relative to the 
P states. The charm quark is too heavy to af- 
fect the running of a at the Bohr radius of T, 
so the spin-averaged spectrum should agree with 
experiment when extrapolated to n/=3. 

The NRQCD and Fermilab groups have also 
calculated spin splittings in the hh system; the 
results are shown in Fig. The (non-quenching) 
systematic errors are expected to be about 5MeV 
for both groups, and are indicated in the figure. 

The most important feature of Fig. Q is that 
the overall size of the Xb splittings agrees with 
experiment, within the expected errors. This is 
significant for two reasons. First, both the Fer- 
milab and NRQCD actions are tadpole improved, 
which increases the strength of the spin interac- 
tions by a factor oiu^^] without tadpole improve- 
ment the Xb sphttings are reduced by 50%[|o|. 
Second, spin splittings are quadratically sensitive 
to the lattice spacing. This is because they are 



roughly inversely proportional to the quark mass. 
Using an incorrectly low value of a^^ to com- 
pare with experiment will result in tuning to too 
heavy a quark mass; this reduces the dimension- 
less spin splittings, which are then multiplied by 
an additional factor of (the incorrect) to make 
them dimensionful. If the correct lattice spacing 
for T physics at /3=6.0 were 2.0GeV rather than 
the 2.4GeV found by the NRQCD collaboration, 
the Xb sphttings in Fig. would be roughly 30% 
smaller in magnitude, resulting in a 15-20MeV 
disagreement with experiment. 

The quenched ratio {xb2 - Xfci)/ (Xfci - X&o) ap- 
pears too small compared to experiment; adding 
dynamical fermions pushes the ratio the wrong 
way. Although this is only a one to two a effect, 
the same behavior occurs with good statistical 
significance in the cc system. In both bb and cc 
the disagreement with experiment is within the 
expected systematic errors. It will be interesting 
to see if improving the NRQCD or IHW actions to 
eliminate 0{v^) corrections to the spin-splittings 
will fix this discrepancy. 



9 



Unfortunately, the ?7f,(lS') has not been ob- 
served experimentally; extrapolating the NRQCD 
results to nf—3 (where a^^ and toq have been 
adjusted) results in a hyperfine splitting of 
44MeV[p^. Because the hyperfine splitting is a 
short-distance quantity, dynamical charm loops 
are likely to have an effect, but a simple pertur- 
bative estimate indicates that this is likely to be 
at most a few MeV. The result should be some- 
where between 45 and 50MeV, which agrees with 
what is expected from potential models. 

5.2. cc 

Both the NRQCD [imi and Fermilab §, 
p6| groups have increased their statistics for 
quenched cc spectroscopy. I do not show the 
spin averaged spectrum, since the only indepen- 
dent level is the 2S, which is too noisy to resolve 
quenching effects. Instead, the spin-dependent 
spectrum is shown in Fig. ||, along with the ex- 
pected (non-quenching) systematic errors (esti- 
mated to be roughly 35MeV), which are much 
larger than for bb since v'^ is three times as big. 
The interesting thing about Fig. || is that these 
systematic errors are dominant; the statistical er- 
rors are negligible. Within the (rather large) sys- 
tematic uncertainty, Fig. || shows agreement be- 
tween simulation and experiment. In addition, 
the trends seen in bb spin splittings (Fig. |^) are 
now quite clear; hyperfine and fine splittings are 
too small and the (xc2 - Xci)/ (Xci - Xco) ratio is 
too big. 

The bad news from Fig. |^ is that cc, for which 
there is more experimental data, is not as well 
controlled a system as bb is. This was to be ex- 
pected; in addition to suffering from worse rel- 
ativistic corrections, cc is a larger meson with 
softer internal glue, and so is more susceptible 
to finite volume effects and the particulars of the 
dynamical fermion content. For these reasons, I 
suspect that the bb system is probably more reli- 
able for as determinations than cc. 

The good news is that cc is a wonderful labo- 
ratory for studying improved actions; the effects 
of (currently neglected) 0{mv^) and 0{asmv'^) 
terms in the NRQCD action are resolvable. A 
drop in the disagreement with experiment of the 
cc spin spectrum of a factor of 3 or so (i.e. w^) 
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Figure 8. cc spin splittings. Symbols are same as 
Fig. ^. Levels are adjusted so that spin averaged 
energies agree with experiment. 



upon inclusion of these terms (or similar ones in 
the IHW approach) would be a tremendous result 
showing the efficacy of improved actions. 

A final comment: the validity of the IHW 
formalism at any value of am is a definite ad- 
vantage for charm physics. As presently formu- 
lated, NRQCD blows up for charm at roughly 
(quenched) /3=5.85, but has a much more easily 
computed improved action. A IHW action, im- 
proved to the some order in and , will almost 
certainly have smaller systematic errors when run 
at a high /3 than the similarly improved NRQCD 
action run at /3maa:=5.85 will have. On the other 
hand, the IHW run (say at /3=6.4) will proba- 
bly cost 100 times as much computer time at the 
same physical volume. Today, with IHW actions 
in use which are less highly improved than the 
NRQCD actions, the relative size of systematic 
errors is unclear. For bb physics, this advantage 
of IHW does not yet enter; even at f3—6A, amwl 
and NRQCD can be used. 
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6. 

Last year, the NRQCD collaboration presented 
a preliminary result for the pole mass of the bot- 
tom quark|3^. The final results of this calcu- 
lation can be found in [ pT[ , which also includes 
the effects of dynamical fermions and a determi- 
nation of the MS mass. This year the Fermilab 
group has presented a preliminary determination 
of the pole mass of the charm quark j|] . 

The basic idea is very simple; given a partic- 
ular regularization scheme (i.e. MS, IHW, or 
NRQCD), the perturbative pole mass can be cal- 
culated in terms of the mass parameter of the 
theory. 



M = Ziml = {\ + C''as{q*))ml 



(9) 



where S labels the scheme, is calculated in 
perturbation theory, and q* is a scale chosen as 
m H or H. The MS mass can be obtained 
by equating the pole mass in the MS and lattice 
schemes; 

^ ' 7MS ^ 7MS ^ ^ ' 

The bare mass of the lattice theory (IHW or 
NRQCD) is determined by requiring that the ki- 
netic mass of the ground state match experiment; 
the result is substituted int o (|9|) to obtain the pole 
mass or ( p^ to obtain the MS mass. 

The NRQCD collaboration obtains a pole mass 
for bottom of Afft=5.0(2) and an MS mass of 

Mj'^S(Af6)=4.O(l)0; the Fermilab group has a 
preliminary result of Mc=1.5(2)[P. 

A word of caution: the NRQCD collaboration 
result has changed since last year. This is because 
q* is actually much softer than the preliminary 
estimate used inQ; ay(q*)«.4 for the bottom 
quark at /3=6.0. This soft scale arises because 
tadpole improvement cancels out the contribu- 
tion to (f of much of the hard glue. A soft q* 
also appears in the MS scheme; applying the pro- 
cedure of to the two-loop MS expression for 
the pole mass[^ gives g*«.22Mb. Because of 
these soft g*s, perturbative uncertainties domi- 
nate the NRQCD determination of Mf,, and make 
an NRQCD determination of impossible. Soft 



g*s are less of a problem in the MS mass deter- 
mination; the soft glue cancels in the ratio of Zs, 
leaving a q* closer to the cutoff. It is not clear 
what non-perturbative effects there might be in 
equating the (perturbative) pole masses of two 
different schemes, especially since infra-red glue 
seems to be playing such an important role. The 
Fermilab group has not found a similar soft q* in 
their IHW calculation. 



Last year, El-Khadra gave an excellent re- 
view of Us determinations from onium spec- 
troscopy [Q, so I will not go into detail. At that 
meeting, there were quenched results from the 
NRQCD and Fermilab groups, as well as prelimi- 
nary dynamical results from the KT group, using 
their nf—2 staggered, /?= 5. 7 ensemble. This year, 
the KTjl|, CDHWQ and NRQCD |l|l groups 
have new dynamical results, using KT , HEM- 
CGC, and HEMCGC ensembles (previously dis- 
cussed), respectively. 

The necessary steps in a lattice determination 
oia^^iMz) are: 

• 1. Determine a{q*) in some lattice scheme; 
in this talk, I will use ay_2(3.41a~^) (de- 
fined by the second order expansion of the 
plaquette in ay )|l|l. This and other pos- 
sible choices will be discussed in depth by 
Michael in his talk||l. 

• 2. Determine a^^ by matching some lat- 
tice "measurement" to experiment. When 
using onium spectroscopy, this is usually 
the spin-averaged IP-IS splitting. The 
NRQCD collaboration has also performed 
their analysis using the T(2S) — T(1S) split- 
ting. Both splittings are very insensitive to 
the bare quark mass and to spin-dependent 
discretization errors. 

• 3. Correct a{q*) to a'"-f^(q*), where n/ is 
the number of quark flavors above thresh- 
old at the typical momentum of the quantity 
used to set the scale in step 2. For example, 
the gluons holding cc or bb together have 
a typical momentum transfer of about 500 
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to 1000 MeV ; the correct physical Uf is 3, 
even though is roughly 8GeV at (3=6.0. 

• 4. Convert from a^"f^{q*) to q;^'(^), us- 
ing a perturbative matching formula. This 
step requires that the a defined in step 1 
above be a good expansion parameter, i.e. 
aiat is not satisfactory. 

• 5. Run a^V) to a^{Mz), using the 3- 
loop beta function. 

The steps above separate the various sources 
of error in the calculation; each step is indepen- 
dent of the others. In particular, the choice of 
scheme for a (step 1) and the choice of experi- 
mental observable to set the scale (step 2) are to- 
tally independent decisions; for example the defi- 
nition of a found in p6| could be combined with a 
scale determination from bb spectroscopy. A non- 
perturbative definition of a {av,2 is such a defini- 
tion, ay determined from the plaquette[^ is not) 
can be used in step 1; only in the conversion to 
MS (step 4) is perturbation theory needed. 

The situation for quenched determinations of 
a^.^{AIz) from onium has not changed since last 
year. Basically, all groups get results consistent 
with .112(8) from cc and .112(5) from 66(see Q 
for details) . This is because the dominant error in 
the calculation comes from step 3, the correction 
to physical fermion content using the technique of 
!§. The KT and NRQCD groups have used this 
technique to correct their nf=2 results to 'T-/=3 
as well, obtaining .112(4) (from cc) and .115(3) 
(from bb), respectively. The agreement between 
cc and bb results is non-trivial; the various deter- 
minations disagree before correction is made for 
quenching. 

This year, a new method for using dynami- 
cal results to correcting quenching effects was in- 
troduced inH^. Results from nf=0 and nf=2 
calculations are run to an arbitrary momentum 
scale near and is extrapolated linearly to 
n/=3. This technique was also used by CDHW. 
The two groups obtain results of .115(2) (from 
bb) and .108(6) (from cc), respectively. The large 
disagreement in central values is due to a dis- 
agreement in the scale determination by the two 
groups, which I discuss below. The NRQCD 



collaboration performed separate, independent, 
analyses for the IP — IS and 2S — IS splittings 
of bb; the central values were almost identical. 
The NRQCD collaboration claims that step 4, the 
perturbative conversion from av,2 to a-^, rather 
than step 3, is now the dominant error in the 
calculation; CDHW, who the Wilson action, cite 
both perturbative conversion and scale determi- 
nation errors as dominant. 

Of all the groups, the NRQCD collaboration 
claims the smallest errors. There are many po- 
tential pitfalls in these as calculations; each of 
them must be addressed before we can be sure 
that the uncertainties are correctly understood. 

The first area of concern is in the choice of 
scheme for a. In order for the calculation to be 
reliable, the conversion from a to a^jg must be 
accurate, i.e. a must be a good perturbative pa- 
rameter. In addition, the observable chosen to 
determine a^^ must exhibit asymptotic scaling in 
a; the ratio aE / (aAa), where E is the observable 
and Aq, is the A parameter of the a scheme, must 
have negligible dependence on (3. This is obvious; 
if it were not the case then determinations of a 
at different f3 using the same observable would 
yield different values of a!^{Mz). Michael will 
compare various prescriptions for a in his talk; I 
discuss asymptotic scaling in the next section. 

A second area of concern is that the dynam- 
ical fermion content could be incorrect. The 
HEMCGC configurations used by the NRQCD 
and CDHW groups had fermion bare masses of 
am=. 01 and am=. 025, i.e. about 25 and 60 MeV, 
while the KT ensemble at am=. 01 corresponds to 
about 30 MeV dynamical quarks. The NRQCD 
collaboration has assumed that the high momen- 
tum transfer inside the T (.5-l.OMeV) will cause 
these these non-physical bare mass values to have 
negligible effect; this should be checked explic- 
itly. In cc, where the momentum transfers are 
not as high, CDHW see a slight dependence on 
the sea quark mass; obtaining a~^'s of 1.90(9) and 
1.72(10) at .01 and .025, respectively Q- There 
are also the usual worries about HMD and the 
'square root trick' for ensembles with two stag- 
gered fiavors. It is important that other calcu- 
lations be done, both with two Wilson and four 



12 



staggered flavors. A four staggered flavor calcu- 
lation would also allow a check of the linear ex- 
trapolation of in Uf. 

The third major source of concern, and the 
most troubling, is the discrepancy between lat- 
tice spacings obtained from NRQCD simulations 
of bb and those obtained by other methods. 

Last year, the NRQCD collaboration reported 
that bb spectroscopy gave a^^~2A GeV at 
quenched P—6.0, significantly larger than typical 
light spectroscopy values of 1.8 - 2.1 GeV. This is 
not alarming since the quenched spectrum should 
not agree with the real world; the f3 function runs 
too quickly so from a high energy observable 
(like bb spectroscopy) should be larger than a^^ 
from a low energy observable (like cc or mp). 

What is disturbing is that dynamical results 
seem to exhibit the same range of a~^'s as 
quenched. The HEMCGC collaboration found 
a very good match between light spectroscopy 
on their (3—5.6, nf~2 staggered ensemble and a 
quenched ensemble at /3«5._95[|^; the NRQCD 
collaboration obtained a bb a^^ of about 2.4 
GeV on the dynamical ensemble, the same as 
at quenched (3=6.0. One would naively expect 
a much smaller spread of momenta at nf=2 if 
quenching were the culprit; perhaps 2.2-2.3 GeV 
when a7r^=2.A GeV . 

bb 

My own (very biased) opinion is that the prob- 
lem probably lies in the light spectroscopy cal- 
culations. This is mainly because we know that 
systematic errors are much more difficult to con- 
trol for light quarks than for heavy; O(oAqc'd) 
scaling violations arising from the Wilson action 
are one obvious candidate. NRQCD, on the other 
hand, appears to have systematic errors under 
control. The spectroscopy results discussed above 
showed no big surprises, and the fact that the 
spin splittings appear to be correct is important, 
since they depend quadratically on the ratio of 
the calculated lattice spacing to its true value. 
The final reason for my bias is that the n/ ex- 
trapolation of NRQCD results is consistent with 
experiment. Fig. ^ implies that the IP-IS and 
2S-1S splittings give different at nj^O but 
agree at n/=3. This is the first spectroscopic cal- 
culation I'm aware of in which quenching effects 
are observed and then removed when dynamical 



fermions are turned on. 

There are three additional calculations which I 
think need to be done before we can be confident 
that the lattice determination of al£l-(Mz) from 
onium really is within 2% of the true value. First, 
the NRQCD calculation using bb spectroscopy 
should be repeated at am=. 025; this will explic- 
itly check the dependence on sea quark mass. Sec- 
ond, a similar calculation must be done using cc 
spectroscopy. This will have to be done at lower 
(3 since NRQCD blows up for charm on the HEM- 
CGC ensemble. If cc spectroscopy yields the same 
value of a at physical nf this will be a beautiful 
example of the removal of quenching effects (since 
a^^s from quenched bb and cc differ by about 
15%). Third, bb and cc calculations must also 
be done at n/=2 using the IHW formulation at 
the .same order as the NRQCD calculation. IHW 
spectroscopy using the tadpole-improved SW ac- 
tion||^ tends to give a~^'s about 10% smaller than 
NRQCD, both for bb and cc; it is important to 
know whether this is because the two actions in 
use have different levels of improvement or if it 
really is a disagreement. At the moment, the 
news here is bad; the obtained from cc by 
CDHW (with no clover term) on the HEMCGC 
configurations is about the same as quenched re- 
sults. It will be interesting to see if a clover IHW 
run on these configurations changes this behavior. 
El-Khadra has estimated the various corrections 
here by using a potential model 

It is important to stress that this lattice spacing 
problem is not just of concern in as calculations. 
For example, both heavy and light quarks appear 
in calculations of /s, discussed by Sommer . 
/b V Mb scales like a~^/^; changing a^^ from 1.8 
to 2.4GeV corresponds to a 40% uncertainty. It 
is very important that we resolve this issue. 

8. Asymptotic Scaling 

As discussed above, a particular observable 
must exhibit asymptotic scaling before it is useful 
in an determination. The V scheme is particu- 
larly convenient for studying asymptotic scaling, 
since aAy is easily determined from the plaque- 
tteQ. Another advantage of expressing lattice 
energies in terms of aAy is the ease of compar- 
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ison between groups; it would be very difficult 
to find highly accurate determinations of nip, for 
example, at the various (3 values I discuss below; 
this is not a problem for the plaquette. 
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Figure 9. Asymptotic scaling plot for the 
quenched bb spectrum. Flat dependence (up to 
logs) on {aAvY indicates asymptotic scaling; lin- 
ear dependence indicates O(a^) discretization er- 
rors. Simulation value of 2S-1S splitting has been 
rescaled by the experimental ratio (1P-1S/2S-1S). 
A perturbative correction has been applied to re- 
move the effects of O(a^) errors in the gluonic 
action. Points at (3=5.7 have been spread hori- 
zontally for clarity. 



The dependence of the spin-averaged quenched 
bb splittings on (aAy)^ is shown in Fig. ^. The 
vertical axis is 

AM _ / {aAM)s^m \ ( 452MeV \ 

Ay ^ aKv ) [(AM),, J ' ^ ' 

where {aAM)sim is the lattice value of the split- 
ting being plotted, {AM)exp is the experimental 
value of that splitting, and 452MeV is the exper- 
imental value of the bb IP-IS splitting (assuming 
a hyperfine splitting of 50MeV) . The rescaling by 
experimental values is equivalent to comparing 
dimensionless ratios to experiment; if two split- 
tings have the same value of aAAI, then their 
ratio matches experiment and they will give the 



same value of a^^. I have chosen (aAy)^ as the 
horizontal axis because the leading discretization 
errors in both the FNAL and NRQCD collabora- 
tion results should be O(a^). All values were ad- 
justed to correct for O(a^) errors in the gluonic 
action. All lattices have physical volume much 
larger than bb mesons. 

There are three different sets of splittings in 
Fig. I (IHW and NRQCD IP-IS and NRQCD 2S- 
IS); all exhibit asymptotic scaling within errors. 
The only point that is inconsistent with flat be- 
havior is the IHW point at /3=5.7; both NRQCD 
splittings at this /3 are consistent with 6.0 (albeit 
with large error bars for 2S-1S). This may be 
because the NRQCD collaboration removes the 
dominant 0{a^) errors in their action (C5 and 
Ce); see § for an estimate of the effects of these 
terms on the Fermilab results. 

Although the results in Fig. ^ exhibit asymp- 
totic scaling, they do not agree on a~^\ The 
discrepancy between IP-IS and 2S-1S within 
NRQCD is good news; this is the quenching effect 
discussed in Fig. ||. The disagreement between 
NRQCD and IHW for the IP-IS splitting will be 
bad news if it doesn't disappear when the IHW 
calculation is repeated with all 0(rnu'*) improve- 
ments turned on. If it doesn't go away and if the 
scale in the NRQCD as calculation is 10% too 
high at ?^/=3, then the value of a^^{Mz) in |l6| 
is shifted down by 2%. 

At low /3, there is an important caveat to the 
above discussion. All of the points in Fig. || had a 
perturbative correction applied to remove the ef- 
fects of O(a^) errors in the gluonic ensembles [p^; 
the uncorrected results are shown in Fig. The 
main effect is to move the (3=5.7 IP-IS points 
down significantly; the (3=5.9 Fermilab point is 
also shifted down about two a. The 2S-1S split- 
ting is negligibly affected, as are cc splittings. 

It is obviously dangerous to rely too heavily on 
points where this correction is large when draw- 
ing conclusions. I have chosen to apply it in these 
scaling plots, however, because we know that it 
is there and a similar calculation of the hyperfine 
splitting works welljl^. It is important to obtain 
direct results using an improved gluonic ensem- 
ble; Lepage showed in his talk that the one-loop 
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Figure 10. Same as Fig. except the perturba- 
tive correction has not been apphed. 



improved gluonic action is tremendously success- 
ful at removing O(a^) errors pof. It will be much 
cheaper to check asymptotic scaling in bb by us- 
ing an improved gluonic action at (3—5.7 than by 
trying to get high statistics at (3—&A. 

Results for the IP-IS splitting in cc are pre- 
sented in Fig. |ll|; 2S-1S results are too noisy to 
be of interest. The discrepancy between NG and 
non-NG splittings in the HS formulation is obvi- 
ous; the two results bracket the IHW value. There 
is only one (3 with NRQCD resuhs; NRQCD 
should blow up at {aKyf < .05. The IHW re- 
sults seem to be scaling better than for bb; see ^ 
for a discussion. The high (3 IHW splittings are 
roughly 10% higher than the low (3 NRQCD split- 
ting. I have no idea if this is related to the similar 
behavior in the bb system; Oimv^) terms left out 
of the NRQCD action should contribute at the 
10% level for cc. It would be very interesting to 
go to the next order of improvement for both the 
NRQCD and IHW actions in this system. 

Note that AA//Ay is roughly 15% larger for 
cc than for bb. Much of this is probably due to 
quenching. Systematic errors in the cc splittings 
could also be noticeably contributing, however. 

I combine cc and bb nf=2 results in Fig. p^. 
Two things are striking. First, the NRQCD bb 
splittings are in much closer agreement than in 
Fig. ^. Again, this the same information as Fig. 



Figure 11. Same as Fig. ||, except displays 
quenched cc spectrum. All splittings shown are 
IP-IS, and have been rescaled by experimental 
ratio between bb to cc. Gluonic correction has 
been applied. NRQCD point at f3=5.7 has been 
moved horizontally for clarity. 



and indicates quenching errors are being removed. 
Second, the NRQCD bb results and CDHW cc re- 
sults disagree by roughly the same amount as in 
Figs. ^ and |ri|. This is the reason for the dis- 
agreement between the two groups' central values 
of ai^(Mz). Unfortunately, the KT group has 
much larger statistical errors. Both the CDHW 
and KT groups used the Wilson, rather than SW 
action; it would be interesting to see the CDHW 
run repeated using a clover, asymmetric clover, 
or 0(mz;^) IHW action. 

9. Conclusions 

We are well on the way to having an excellent 
understanding of HH mesons. We already have 
many cross-checks; more should be forthcoming. 

I wish to thank my colleagues in the Fermilab, 
Kyoto-Tsukuba, CDHW and NRQCD collabora- 
tions for their patience and work in supplying me 
with and explaining their results. This talk was 
prepared at the University of Glasgow, during a 
visit supported by the PPARC. I am extremely 
grateful to the members of the Physics Depart- 
ment there, especially Christine Davies, for their 
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Figure 12. Same as Fig. ^, except displays nf=2 
spectrum. NRQCD points are bb, all others are 
cc. Splittings have been rescaled by experimental 
ratio. NRQCD IP-IS and KT staggered points 
have been moved horizontally for clarity. KT 
staggered points are a lower bound, due to finite 
volume effects. 



hospitality during my stay. This work was sup- 
ported in part by a DOE grant. 
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